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Abstract 

Imposing Huygens' Principle in a 4D Wightman QFT puts strong con- 
straints on its algebraic and analytic structure. These are best understood in 
terms of "biharmonic fields" , whose properties reflect the presence of infinitely 
many conserved tensor currents. In particular, a universal third-order partial 
differential equation is derived for the most singular parts of connected scalar 
correlation functions. This PDE gives rise to novel restrictions on the pole 
structure of higher correlation functions. An example of a six-point function is 
presented that cannot arise from free fields. This example is exploited to study 
the locality properties of biharmonic fields. 

1 Introduction 

In spite of the quantitative successes of renormalized perturbation theory, no non- 
trivial quantum field theory (QFT) in four spacetime dimensions (4D) has been 
constructed rigorously. It is sometimes suggested that the Wightman axioms might 
be too restrictive, but most attempts at relaxing them lead to physically unac- 
ceptable consequences. E.g., a violation of causality cannot be exponentially small 
|14j . Admitting indefinite (physical) Hilbert spaces not only jeopardizes the statis- 
tical interpretation of correlations, but it makes statements about convergence of 
approximations delicate if not meaningless. 



'Lecture given by K.-H.R. at the Workshop "Lie Theory and Its Applications in Physics", 18-24 
June 2007, Varna, Bulgaria; to appear in the proceedings 
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On the contrary, the usual attitude is to strengthen the axioms, e.g., by imposing 
additional symmetries or phase space or additivity properties. This facilitates the 
analysis of models but clearly reduces the number of theories. Even worse, demand- 
ing too much may trivialize the theory. E.g., requiring that the two-point function 
be supported on a mass shell [9] or only to have finite mass support [5], already 
forces the field to be a (generalized) free field. 

It is one of the merits of the axiomatic approach [18, 8j that it allows to pinpoint 
such obstructions even before a model is formulated. It makes no assumptions about 
intermediate steps and limits through which a theory is constructed. Referring only 
to its intrinsic features, it avoids assigning significance to the artifacts of the descrip- 
tion. Especially, it guides navigation between Scylla (physical meaningless theories) 
and Charybdis (free fields) in the quest for a proper set of model assumptions that 
can possibly be satisfied by other than free fields. 

In a series of recent papers [121 US], two of us have suggested to demand "global 
conformal invariance" (GCI) which is the postulate that the conformal group is rep- 
resented in a true (i.e., not a covering) representation. (The term "global" stresses 
the fact that the conformal transformations are defined globally on the Dirac com- 
pactification of Minkowski spacetime.) 

GCI implies rationality of all correlation functions |13j . On one side, this is 
a desired feature since it allows to parameterize each correlation function by the 
coefficients of a finite set of admitted rational structures. The latter are determined 
by conformal invariance, and by the unitarity bound on the representations of the 
conformal group that possibly contribute to operator product expansions (OPE) 
giving rise to upper bounds on the pole in each pair of variables. The coefficients 
are then further restricted only by Hilbert space positivity. While the latter is 
highly nontrivial to control in general, there has been considerable progress at the 
four-point level [10] . 

On the other side, GCI is equivalent to Huygens locality |10j . i.e., commutativity 
of fields not only at spacelike but also at timelike separation. This feature seems 
to be conspicuosly close to free field theory, since any scattering of field quanta 
should give rise to causal propagation within the forward lightcone. Indeed, it has 
been shown [4] that if a Huygens local QFT has a complete particle interpretation, 
then its scattering matrix is trivial. Therefore, GCI allows only nontrivial theories 
without asymptotic completeness, which may not be entirely unphysical in a scale 
invariant theory. Note that by rationality, all scaling dimensions are integer numbers, 
which does not mean that they are canonical. It is conceivable that in a QFT 
with anomalous dimensions, at some value of the coupling constant all dimensions 
simultaneously become integers. Examples for such a hazard are well known in 2D 
conformal QFT [6]. 

We shall report here on a recent analysis [TT] of the consequences of the following 
fact. For a symmetric tensor field of rank r and dimension d, one calls d — r the 
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"twist". Twist two fields are necessarily conserved, because their two-point func- 
tions, completely determined by r and d, are conserved. (Hilbert space positivity 
crucially enters in this argument: if the norm square of a vector vanishes, then the 
vector itself vanishes. The rest of the argument invokes the Reeh-Schlieder theorem 
for which in turn locality and energy positivity are essential.) These conservation 
laws can be reformulated in the form 

D x V(x,y)=0 = D y V(x,y) (1.1) 

where V{x,y) is the sum of all twist two contributions to the OPE <f>i(x)(f>2(y) of 
two scalar fields of equal dimension d. We say that the "bi-field" V is biharmonic. 

Our first main result is the unravelling of a hidden consequence of this equation: 
a third-order partial differential equation to be satisfied by (certain parts of) all cor- 
relation functions involving <f>i{x)(p2{y)i that is a necessary and sufficient condition 
for biharmonicity. 

We present a solution to this PDE, and the corresponding transcendental six- 
point correlation function of V(x, y) that cannot be produced by Wick products of 
free fields. On the basis of this solution, which we believe to be prototypical for the 
general case, we then study the locality properties of the bi-field. 

2 Biharmonicity 

We first explain how harmonicity of V(x, y) serves to define its correlation functions. 
Because V(x, y) is biharmonic, there are in fact two such prescriptions, and hence 
V(x, y) is overdetermined. The consistency condition gives rise to a restriction on 
the correlation functions involving <j)\{x)4>2{y)- 

The crucial fact is the following "classical" result [3]. If u(x) is a power series in 
x G R or C D , then there is a unique "harmonic decomposition" 

u(x) = v(x) + x 2 u(x), (2-1) 

such that v is harmonic (O x v = 0) and u is again a power series. We call v the 
"harmonic part" of u. (Questions of convergence will be discussed in Sect. 5.) 

We apply this fact as follows. Twist two contributions to correlations involving 
<f)\(x)(j)2{y) have a leading singularity ((x — y) 2 ) 1 ~ d , while all higher twist contri- 
butions are less singular. Therefore, the (Huygens bilocal, but not conformally 
covariant) bi-field U (x, y) defined by 

M*)Mv) ~ (Mx)Mv)) =■ ((* - yfr {d ~ 1] ■ u(x, y ) (2.2) 

is regular at {x — y) 2 = 0. U (x, y) contains all twist > two contributions to the 
OPE, and because the twist > two contributions are suppressed by another factor 
{x — y) 2 , we may write 

U(x,y) =V(x,y) + (x-y) 2 U(x,y) (2.3) 
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where both V and U are regular at (x — y) 2 = 0. Now consider any correlation 
function 

u(x, y, . . . ) = (U(x, y)4> 3 (x 3 ) ■ ■ ■ (p n (x n )), (2.4) 

where . . . stands for the arguments of all other fields. Its Taylor expansion in x 
around y is a power series in x — y with coefficients independent of x. Thus, by 
(|2.3p and because \I\ x V(x,y) = 0, the desired correlation function (again as a power 
series) 

v(x, y, . . . ) = (V(x, y)<t>z{xz) ■ • • 4>n{x n )} (2.5) 

is the harmonic part of this series. By construction, v (x, y, . . . ) transforms like the 
correlation function of a conformal bi-scalar of dimension (1, 1). 

On the other hand, the Taylor expansion of u(x, y, . . . ) in y around x is another 
power series in x — y, whose coefficients do not depend on y, and because also 
\Z\ y V(x, y) = 0, v(x, y, . . . ) may as well be determined as the harmonic part of this 
latter series. 

This overdetermination imposes a consistency condition on the function u(x, y, . . .). 
Its nontriviality may be seen from the following example. Consider u(x,y, . . . ) = 
(y — xq) 2 /(x — x 3 ) 2 {y — Xi) 2 (y — x^) 2 . This function is harmonic with respect to 
x, hence is harmonic part with the first prescription coincides with u itself, but it 
is not harmonic with respect to y, so the harmonic part with respect to the second 
prescription differs from u, and the two definitions of v are conflicting each other. 
Thus, a function u(x, y, . . . ) as in this example cannot occur as a correlation function 
of U(x, y). 

We conclude that biharmonicity of the bi-field V(x,y), which follows from the 
conservation of conformal twist two tensor fields, implies a nontrivial restriction on 
the possible correlation functions u(x,y,...) of the bi-field U(x,y), and hence on 
the correlations involving (pi(x)<j)2(y)- 

We shall now turn this condition into a partial differential equation. 

Global conformal invariance implies that scalar correlation functions are Laurent 
polynomials in the variables = pji = (xi — xj) 2 of the form 

(Mxi)---Mxn)) = Yl f J u ( 2 - 6 ) 

where the integer powers fiij = fiji satisfy the homogeneity rules 

^2.fJ<ij = -di, (2.7) 

and the absence of non-unitary representations of the conformal group in the OPE 
implies the lower bound for the connected parts of (|2.6[) 



2^/^' ~~^ > d<i dj ' 



(2.8) 
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Let 0i and 4>2 have the same dimension d. It follows that all correlations (|2.4p (that 
give contributions to (|2.6p ) are Laurent polynomials in pij, which are separately 
homogeneous of total degree —1 in pi^ {k ^ 1) and in p2k (k ^ 2), and which are 
true polynomials in pi 2 . Because all terms involving a factor p\2 have zero harmonic 
part in the harmonic decompositions, we need to consider only the function uq which 
is the contribution of order (P12) of u. Then uq is separately homogeneous of total 
degree —1 in pu- > 2) and in p2k (k > 2). 

It is now important that the harmonic part v is a real analytic function in a 
neighborhood of x\ = x 2 , provided (x2 — X{) 2 / for all j > 2 (see Sect. 5). We 
may therefore expand v as a power series Y^=i h n /n\ ■ p\ 2 . The coefficients h n are 
functions of all the remaining variables p™ ^ p±2, and are separately homogeneous 
of total degree — n — 1 in p\^ (k > 2) and in p 2 ft (k > 2). 

Let us write djk = dkj = Then the wave operator H Xl has the form 

D -i = - A {T,2< j<k <n P ^ dlk ) = " 4 Pl +^1^12) (2-9) 

valid on homogeneous functions of total degree —1 in pi^ (k ^ 1) [TO], where 

D i = 'Y]„^. , ^ Pjkdijdik and E 1 = VV. P2idu- (2.10) 

Similarly, replacing 1 <-> 2 everywhere, one represents = — 4(L>2 + -E-2 $12)- Thus, 
the two conditions n xi / = = C\ X2 f give rise to two recursive systems of partial 
differential equations for the coefficient functions h n of the form 

Ei h n+ i = -Di h n and E 2 h n+1 = -D 2 h n . ( 2 -H) 

At n = we obtain the integrability condition 

{EkD 2 - E 2 D X ) h = (E 2 E 1 - E X E 2 ) h x . (2.12) 

Because h\ is separately homogeneous of total degree —2 in pik (k > 3) and in p 2 k 
(k > 3), the commutator [E 2 E\ — E\E 2 ) vanishes on h±. Since v is the harmonic 
part of u, its leading term h,Q equals uq. Thus, we arrive at 

Result 1: The function uq solves the third-order partial differential equation 

(E 1 D 2 -E 2 D 1 )u = 0. (2.13) 

Next, when (|2.13p holds and the recursion is solved for h\, one has (D\E 2 — 
D 2 E\)h\ = —{DiD 2 — D 2 Di)ho = because D\ and D2 commute. But D\E 2 — 
D 2 E\ = E2D1 — E\D 2 . Therefore, the integrability condition for the next step of 
the recursion is automatically satisfied, and the argument passes to all the higher 
steps. Thus, (|2.13p secures solvability of the entire recursive systems (|2.1ip . 
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3 Consequences 

One should worry how there can be a new universal (model-independent) partial 
differential equation for the correlation functions. It is important to notice that this 
PDE can not be regarded as some "equation of motion". The reason is that it is 
satisfied only by the leading part uq of n. The splitting of u into uq+ a remainder 
does not correspond to any local decomposition of the bi-field U(x,y). Thus, the 
PDE (|2.13p cannot be formulated as a differential equation for some (bi-)fields in 
the theory. 

Instead, it should be understood as a kinematical constraint. Because its solu- 
tions no must at the same time be Laurent polynomials, the PDE rather selects a 
(finite) set of admissible singularity structures, that depends on the dimensions of 
the scalar fields involved through the lower bounds on pij. 

Indeed, we have shown in [TT] that the PDE (j2. 13[) implies the following constraint 
on the pole structure of a Laurent polynomial no in pn- and p2k (k > 2), that is 
homogeneous of degree —1 in both sets of variables separately: Suppose no contains 
a monomial 

IL >2 PlkP2k X other factors I 3 - 1 ) 

where the other factors depend only on pj~i (k, I > 2). If there are i ^ j such that 
pu < and < (a "double pole" in xi), then one must have p2k > for all 
k > 2, k 7^ In particular, this excludes "triple poles", because a triple pole 
in x\ would imply that all p2k > 0, contradicting homogeneity. The most involved 
possible pole structure of no is therefore of the form 

polynomial 

— p — j — - — — x other factors (3.2) 
PliPljPliPlj 

where the polynomial takes care of the proper homogeneity. The corresponding 
contribution to the connected correlations involving (pi(x)4>2(y) is therefore 

1 polynomial , . 

x other factors. (3-<j) 



d—l ,3? „1 „r „s 
Pl2 PxiPxjPliPlj 

Note that no constraints arise on higher twist contributions {p\2 > 1 — d) or on poles 
of two fields with d\ ^ d.2- 

The interest in double poles is due to the fact that twist two bi-fields made of free 
fields, such as :(p(x)(p(y): or [x — ?/) A ':'0(x)7 /i V ; (y) : are always Wick bilinears, so that 
their correlation functions can never contain a double pole. A nontrivial double pole 
solution (an example will be displayed below) is therefore a candidate for a Huygens 
local QFT not generated by free fields. 

We have also shown in [TT] that 
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Result 2: A correlation function involving V{x\,X2), i.e., the harmonic part of the 
Laurent polynomial uq{x\,X2, ■ ■ ■), is again a Laurent polynomial if and only if uq 
does not contain a double pole in x\ or X2- 

Four-point functions (U(xi, #2 )</>3 (£3)04(^4)) can never exhibit double poles in 
X\ or X2, just "by lack of independent variables". Therefore, four-point functions 
of twist two bi-fields are always rational. From this one can deduce that their 
partial wave expansion cannot terminate after finitely many terms, i.e., the OPE of 
4>i(%)4>2(y) must contain infinitely many conserved tensor fields. 

If all fields are scalar s of dimension 2, hence nij > — 1, double poles cannot occur 
in any n-point function subject to the cluster decay property. We have exploited 
this fact in [11] to prove that scalar fields <j) of dimension 2 are always Wick products 
of the form Y^-^ij :L Pi L Pj'-{ x ) °f massless free fields. In this argument, Hilbert space 
positivity plays a crucial role because one has to solve a moment problem in order 
to get the correct coefficients for all n-point functions simultaneously. (When we 
do not insist that the theory possesses a stress-energy tensor with a finite two-point 
function, then the fields (f> may also have contributions of generalized free fields.) 

The simple pole structure of correlation functions of dimension 2 fields can be 
converted into commutation relations of the twist two biharmonic fields occurring in 
their OPE. The result is an infinite-dimensional Lie algebra, whose unitary positive- 
energy representations can be studied with methods of highest weight modules. It 
turns out that there are no other representations than those induced by the free field 
construction [2]. 

4 An example with double poles 

The following six-point structure solves the PDE (|2.13|) both in the variables x\, X2 
and in the variables X5 , xq : 

. . (P15P26P34 - 2^i5P23P46 - 2pi5P24P36) [121 [5,61 , , _ s 

«(«!,..., a 6 ) = , (4.1) 

P13P14P23P24 • P34 • P35P45/536P46 

where (. . . stands for the antisymmetrization in the arguments Xi, xj, and pij = 
[x{ — Xj) 2 as before. This structure in addition obeys all homogeneity rules (|2.7p . 
pole bounds (|2.8p and cluster conditions in order to qualify as (a contribution to) 
the correlation function 

(U( Xl , X2W{x^)4>'{x A )U(x 5 , x 6 )) (4.2) 

where the scalar field <j)' has dimension d! . The multiple poles in the variables x$, X4 
do not contradict the previous argument (Sect. 3) excluding triple poles in the twist 
two "channel" , when either d (the dimension of the fields <j)\ , 4>2 in (|2.2p , generating 
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U) or d! is > 2, because they don't arise in a channel of twist two (I/P34 3 is twist 
six, and l/psi and l/p4t are twist two only if d = d! = 2). 

We determined the corresponding (contribution to the) correlation 

(V(xi, x 2 )4>'{x 3 )4>'(x i )V(x 5 , x 6 )), (4.3) 
• • • , Xq), as the (simultaneous) harmonic part(s) of u{x\, . . . , xq). Let 



P12P34 

Pl3/>24 



t 



P14P23 
P13P24 



(4.4) 



denote the conformal cross ratios, and s' and t' the same with 1, 2 replaced by 5, 6. 
Then 



v(xi, . . . , x 6 ) = u(xi, . . . , x e ) ■ g{t, s)g(t', s') + 

2 (P13P24 • P35P46)[i,2][5,6] 



d'— 2 

P13P14P23P24 • P34 • P35P45P36P46 



l- ff (t,s) 5 (f,/)) 



(4.5) 



has the required power series expansion u{x\, . . . ,xq) + 0(/?i2>P56) provided g(s,t) 
is of the form g(s,t) = J2 n >o 9n(t)/nl ■ s n with go(t) = 1, and it is harmonic in all 
four variables x±, x 2l £5, xq provided g solves the PDE 

((1 - td t )(i + td t + sd s ) - [(1 - td t ) + 1(2 + tdt + sd s )]d s ) g = o. (4.6) 

The solution is 



+ 



Li 2 (u) + Li 2 (v) — Li 2 {u + v — uv) 
-u \ _ . / 



+ 



(4.7) 



Li 2 



1 - u 



1-u 



Li 2 



uv — u — V 
(l-u)(l-v) 



where it and u (apologies for the duplicate use of letters!) here stand for the "crural" 
variables defined by the algebraic equations 



uv 



and 



t = (1 -«)(!-«). 



(4.8) 



Li 2 is the dilogarithmic function defined by analytic continuation of its integral or 
power series representations (0 < x < 1) 



Li 2 (x) 



Mi - 1) 



n>0 



(4.9) 



Notice that g is regular at s = in spite of the prefactors ~ 1/s. This transcendental 
correlation function can definitely not be produced by free fields. It was found by 
turning the differential equation (|4.6p into the recursive system 



(1 + (n + l)t - t(l - t)d t )g n = (1 - td t ){n + td t )g n -i 



(4.10) 
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with go(t) = 1, and resumming the solution 




■ 2 F 1 {n,n + l;2n + 2;l-t) 



(4.11) 



by exploiting the integral representation of hypergeometric functions. 

5 Local commutativity 

We shall now discuss the issue of local commutativity of the bi-field V(x,y). The 
naive argument would go as follows: since U(x, y) is Huygens bilocal in the sense of 
local commutativity for spacelike or timelike separation from x and y, the correlation 
functions 



are independent of the position k where U(x, y) is inserted. By the uniqueness of 
the harmonic decomposition, the same should be true for their harmonic parts 



hence V(x,y) commutes with (f> k (xk)- 

However, this argument is not correct because of convergence problems of the 
power series. The transcendentality of the correlation function (|4.5j) shows that 
V(x, y) in this case is certainly not a Huygens bilocal field, which must have rational 
correlation functions by the same argument [13] as for Huygens local fields. On 
the other hand, the Result 2 in Sect. 3 yields a necessary and sufficient condition 
(obviously violated by (|4.1|) ): 

Result 3: V(x, y) is Huygens bilocal if and only if the coefficients of the twist two 
pole {(x — y) 2 )^^ 1 ^ in every correlation involving <p(x)4>(y) (i.e., the leading parts 
uo of the correlations ofU(x,y)) never exhibit "double poles" in the variables x or 
y (as explained Sect. 3). 

In general, i.e., when there are double poles, V(x,y) is originally only defined as 
a formal power series (in x — y) within each correlation function. Even when these 
series converge, it is not a priori clear what the labelling pair of points x, y has to do 
with its localization in the sense of local commutativity with other fields, because 
splitting the OPE into pieces is a highly nonlocal operation (involving projections 
onto eigenspaces of conformal Casimir operators). 

In order to study local commutativity, we need to control convergence of the 
series defining the harmonic part. The latter can be addressed with the help of the 
"generalized residue formula". This integral representation of the harmonic part 
was found recently [T] in the context of higher-dimensional vertex algebras: 



u k (x,y, ■■■) = (<p3{x 3 ) . . . 4>k (x k )U(x,y)<j> k +i (x k+ 1 ) . . . 4> n (x n )) 



(5.1) 



v k (x, y,...) = (03 (a*) • • • (f>k{xk)V(x, y)4> k +i(xk+i) ■ ■ ■ (t>n(x n )) 



(5.2) 




(5.3) 
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Here, z 2 = X^£=i z a * s the complex Euclidean square. M r is the compact submanifold 
r • S 1 • S^ -1 (S 1 C C is the complex unit circle, and S- 0-1 C M° C C D the real unit 
sphere), and d D z\m t the induced complex measure. 

The radius r > has to be chosen such that u(z) converges absolutely for z £ M r . 
Then for x small enough such that the kernel converges absolutely as a power series in 
x for every z £ M r , the integral converges as a power series in x, and is independent 
of the choice of r. 

This formula for the harmonic part w.r.t. the Euclidean Laplacian remains valid 
for the Lorentzian Laplacian, provided z 2 is replaced by the (complex) Lorentzian 
square, and the unit sphere by the set {(ix°,x) : (x°,x G S D_1 )}. This is true 
because the map C D — ► C D , (z°,z) i— > (iz°,z), intertwines the Euclidean with the 
Lorentzian harmonic decomposition. 

In the case at hand, where x — y plays the role of x and u(x) is the Taylor 
series around x or around y, respectively, of a Laurent polynomial with poles at 
(x — Xj) 2 = and (y — Xj) 2 = 0, we find absolute convergence in the domain 



in the first case, and the same with x replaced by y in the second case. Especially, if 
x and y are spacelike or timelike separated from all other points Xj, these domains 
are not empty. We have therefore 

Result 4: The formal power series Vk(x,y, . . . ) for the correlation functions ij5.2j) 
converge absolutely within the domains i5.4\) , and the resulting functions Vk do not 
depend on the position k where V(x, y) is inserted in i5.2]) . 

The issue of local commutativity of V(x,y) with <frk(xk) now amounts to the 
question whether Vk-i = still holds outside the domain (|5.4p . as long as x^ is 
spacelike (or timelike) from x and y. We conservatively anticipate that the correla- 
tion functions are real analytic functions of real spacetime points within the region 
where local commutativity holds. Then the existence of a unique real analytic con- 
tinuation from (|5.4p to some other configuration implies Vk-i = Vk at the latter 
configuration by virtue of Result 4, and hence commutativity. Continuation beyond 
a singularity requires to go through a suitable complex cone which depends on the 
position k where V(x,y) is inserted in (|5.2p . hence commutativity will fail. Put 
differently, our strategy to establish locality by inspection of analyticity inverts the 
usual axiomatic reasoning [18] by which one derives the domain of analyticity from 
the known locality (and energy positivity). 

We want to discuss specifically the local commutativity of V(xi, 22) with <j)'(xz) 
in the case of the example (14. 5p . by studying its maximal real analytic continuation 
starting from the domain (|5.4p . which is a neighborhood of x\ = X2 where s = 0, t = 
1, hence u = v = 0. Clearly, we can only reach configurations where (x\ — Xk) 2 7^ 



x-y\\ + y/\\x - y\\ 2 + \(x - y) 2 \ < 




n 



(5.4) 
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has the same sign as (x2 — x^) 2 for k = 3 and k = 4, because this is trivially true at 
x\ = X2 and we cannot pass through t = or t = oo where the variables u or v in 
(|4.7p would hit the singularities of the dilogarithmic function Li2(z) at z = 1 and 
z = oo. 

We claim that (|4.5j) has a unique real analytic continuation to all these points, 
or equivalently, that g(s,t) given by (|4.7p has a unique real analytic continuation 
in the region t > 0, s arbitrary (real). This is obvious for the last terms in the 
two lines of (|4.7p because for t > their arguments are < 1. For the study of the 
remaining terms, we solve (|4.8[) for it and v (where it does not matter which one is 
which because of the manifest symmetry of (|4.7p under u *-* v) 

u, v = - (l - t + s ± V(l - t + s) 2 - 4s) . (5.5) 

In the range s < (1 - \/i) 2 , u and ■y are real and u + v = 1 — t + s < 2. From 
(l—u)(l—v) = t > 0, we see that both u and t> are < 1, and so are an d i"r^- The 
continuation to these points is unambiguous. In the range (1 — Vt) 2 < s < (l + ^/i) 2 , 
u and v are complex and conjugate to each other, so that the first two terms in both 
lines of (|4.7p are always the sum of the values on the two branches above and below 
the cut. In particular, g(s,t) is real and L12 in (14. 7p may be replaced by its real 
part. Finally, in the range s > (1 + Vi) 2 , we find u + v > 2, hence both u and v and 
a ls° j^z and are > 1. All four arguments hit the cut of L%2(z). But because its 
discontinuity is imaginary, the real parts are real analytic. This proves the claim. 

As explained before, the maximal domain of real analyticity specifies those con- 
figurations xi,X2,x$, where <f>'{x$) commutes with V(x\,X2)- We may assume 
x\ — ► ±00 (which can be achieved by a conformal transformation), hence t = 
(X2 — X3) 2 /(xi — X3) 2 . Thus we get commutativity whenever X3 is simultaneously 
spacelike or timelike from x\ and X2- We summarize 

Result 5: The transcendental (part of a) correlation function &4.5\) is compatible 
with local commutativity between V(x, y) and 4>'(z) when x — z and y — z are either 
both spacelike or both timelike. 

The set of these configurations is locally, but not globally conformal invariant, 
since a conformal transformation may switch the sign a of (x — z) 2 /(y — z) 2 . This is 
not a contradiction: connecting configurations with a = + with those with a = — by 
a path in the conformal group, one must necessarily pass through x = 00 or y = 00, 
where the OPE in terms of power series ceases to make sense. This breakdown of 
GCI for the biharmonic field V(x,y) is, of course, just another manifestation of its 
violation of Huygens bilocality. 

It is worth noticing that another decomposition theory for the OPE in conformal 
QFT was developed in [T6], [T7] . While it is coarser than the twist decomposition (it 
is even trivial in the GCI case), it was found to exhibit, at least in two dimensions 
|15j . a "localization between the points" with similar implications as the present 
one. 
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6 Conclusion 

We have outlined recent progress in the intrinsic structure analysis of quantum field 
theories in four dimensions, under the assumption of "global conformal invariance" 

We have found nontrivial restrictions on the singularity structure of correlation 
functions. The encouraging aspect is that these restrictions allow a small "margin" 
beyond free correlations, for which we have given a nontrivial example. It exhibits 
a local but not Huygens local bi-field V{x, y) whose correlation functions involve 
dilogarithmic functions. Local commutativity with a third field at a point z is 
shown to hold (in this example) whenever x — z and y — z are either both spacelike 
or both timelike. The possible failure of local commutativity when one is spacelike 
while the other is timelike, occurs only in correlations of at least five points, because 
four-point functions cannot exhibit the characteristic "double poles" in the twist two 
channel that are responsible for the transcendental correlations involving V(x,y). 

A serious question remains to be answered before our six-point structure is estab- 
lished as (a contribution to) a manifestly non-free correlation function: we cannot 
control (at the moment) Wightman positivity at the six-point level. In the case at 
hand, this means that we do not know whether the vectors (p'(z)V(x, y)\0) span a 
Hilbert (sub-)space with positive metric. Because our six-point structure reduces 
in the leading OPE channels to five- and four-point structures that can also be 
obtained from free fields, positivity can only be violated in higher channels where 
our present knowledge of partial waves is not sufficient. Far more ambitious is the 
problem whether a given six-point function can be supplemented by higher correla- 
tions satisfying Wightman positivity in full generality, i.e., to recover the full Hilbert 
space on which the bi-field V(x,y) and its generating field 4>{x) act. 
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